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Recently Magnea [Phys. Rev. D61, 056005 (2000); Phys. Rev. D62, 016005 (2000)] claimed to 
have computed the first sum rules for Dirac operators in 3D gauge theories from OD non-linear o 
models. 1 point out that these computations are incorrect, and that they contradict with the exact 
results for the spectral densities unambiguously derived from random matrix theory by Nagao and 
myself. 
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Magnea has recently claimed to have derived Dirac 
spectral sum rules for three-dimensional gauge theories 
coupled in a (P, Z2)-invariant manner to fundamental 
fermions with Nc — 2 (corresponding to the Dyson in- 
dex /3 = 1) and adjoint fermions (/3 = 4). She employed 
the small-mass expansion of the low-energy effective the- 
ories, i.e. the zero-dimensional tr models over Rieman- 
nian symmetric spaces M. = CII and BDI, instead of 
AlII that had been proposed for the case with funda- 
mental fermions and JVc > 3 (/3 = 2) ||^. She concluded 
that the first sum rule in the presence of even number 
(Nf) of massless 2-component complex or 4-component 
real (Majorana) fermions is common both to the (3 — 1 
and /? = 4 universality classes, and takes the form^^ 
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where ( stands for an unfolded Dirac eigenvalue (i.e. 
rescaled by l/{Trp{0))) and pif^\() stands for the scaled 
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spectral density. If true, this conclusion, derived from an 
obviously correct formula (see Ref. |3|) 
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(d stands for the rank of the matrix that parameterizes 
M, and M for the dimension of A^), would be surprising, 
as the four-dimensional counterpart of the spectral sum 
rules are known to be different for three values of /3 . 

On the other hand, Nagao and myself have ob- 
tained Pfaffian expressions for the generic p-level correla- 
tion functions in a presence of an arbitrary number of (for 
f3 — 1) and an arbitrary number of pairwise degenerate 
(for (3 — 4:) finite fcrmion masses {/i/}, by applying the 
skew-orthogonal polynomial method to pertinent random 
matrix ensembles. To make comparison with Eq.([^), I 
take a completely confiuent limit /iy — > for all /, of our 
results [Ref. g, Eqs.(2.40), (2.42), (3.23), (3.25), (3.49), 
(3.51)] with p — I (spectral density), to obtain: 
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(Si stands for the sine-integral function) and so forth. 
These expressions for the spectral densities lead to the 
sum rules 
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which are sensitive to the Dyson index /3. The above sum 
rules agree perfectly with the numerical results for ran- 
dom matrix ensembles of large but finite ranks (~ 40), 
obtained by Hilmoine and Niclasen via two alterna- 
tive methods (an analytical method of Widom's [Table 
4 of Ref. 1^] and numerical Monte-Carlo simulations of 
random matrix ensembles). Therefore I conclude that 
the expression (|l|) is erroneous, and the coincidence of 
the first sum rules for (3 = 1 and /3 = 4 claimed in her 
papers is illusory. 
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Nf in Ref. ||^ stands for the number of 4-component 
spinors, i.e. a half of Nf in this note and in Refs. gj^,^. 
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